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The existence and nature of tripartite entanglement of a noninteracting Fermi gas (NIFG) is 
investigated. Three new classes of parameterized entanglement witnesses (EWs) are introduced 
with the aim of detecting genuine tripartite entanglement in the three-body reduced density matrix 
and discriminating between the presence of the two types of genuine tripartite entanglement, W\B 
and GHZ \ W . By choosing appropriate EW operators, the problem of finding GHZ and W EWs 
is reduced to linear programming. Specifically, we devise new W EWs based on a spin-chain model 
with periodic boundary conditions, and we construct a class of parametrized GHZ EWs by linearly 
combining projection operators corresponding to all the different state-vector types arising for a 
three-fermion system. A third class of EWs is provided by a GHZ stabilizer operator capable of 
distinguishing W\B from GHZ\B entanglement, which is not possible with W EWs. Implementing 
these classes of EWs, it is found that all states containing genuine tripartite entanglement are of 
W type, and hence states containing GHZ \ W genuine tripartite entanglement do not arise. Some 
genuine tripartite entangled states that have a positive partial transpose (PPT) with respect to some 
bipartition are detected. Finally, it is demonstrated that a NIFG does not exhibit "pure" W \ B 
genuine tripartite entanglement: three-party entanglement without any separable or biseparable 
admixture does not occur. 

PACS numbers: 03.67.Mn, 03.65.Ud, 71.10.Ca 



I. INTRODUCTION 

The existence of entangled states is a distinctive sig- 
nature of quantum mechanics of the most profound con- 
ceptual and practical importance. The phenomenon of 
entanglement promises to be the source of many far- 
reaching technological advances of the 21st century. A 
fundamental and quantitative understanding of its na- 
ture will be instrumental to its exploitation in quantum 
information processing based on optical and condensed- 
matter systems. Moreover, such an understanding can 
bring new insights into the microscopic physics, statisti- 
cal mechanics, and phenomenology of strongly interact- 
ing quantum many-body systems. 

There is a growing body of work that seeks to estab- 
lish the entanglement properties of the states of quantum 
many-body systems and the roles entanglement plays in 
the observed behavior of these systems. In particular, 
much progress has been made on the entanglement prop- 
erties of spin-lattice models, stimulated by the pioneering 
studies of Osterloh et ai. [l[ and Osborne and Nielsen Q . 
The recent review of Amico et al. [|| assesses the state 
of knowledge on bipartite and multipartite entanglement 
for diverse many-body systems including spin, fermion, 
and boson models. The great majority of examples stud- 
ied involve spin systems and systems of particles made 
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distinguishable by localization. 

There is an ongoing debate on the nature of entangle- 
ment in systems of identical particles — just how does the 
indistinguishability of the particles impact the quantifi- 
cation of entanglement? Attempts at clarification of the 
various issues that arise when different Bose and Fermi 
systems with different degrees of freedom are studied 
|3l— ll 31] has led to the examination of various quantities 
deemed to measure or detect entanglement in the pres- 
ence of indistinguishability. Here we shall (i) focus on 
the noninteracting Fermi gas (NIFG) as represented by 
the three-fermion reduced density matrix of its ground 
state and (ii) adopt the entanglement witness (EW) cri- 
terion [T3, [lf| for analysis of the entanglement content 
of this tripartite state descriptor, which overcomes disad- 
vantages of some of the earlier work. We shall introduce 
new classes of parameterized EW operators for indistin- 
guishable fcrmions, to enable detection of GHZ\W (i.e., 
the subset of GHZ subtracted by the set of W) and W\B 
genuine tripartite entanglement in the NIFG, if one or the 
other is present. 

This paper is organized as follows: In Sec. II we review 
certain definitions basic to the discussion of entanglement 
in systems of distinguishable or indistinguishable parti- 
cles, and some existing formalism and results from other 
authors relevant to characterization of the entanglement 
properties of the noninteracting Fermi gas [1, 0, lii-fllj . 
Importantly, we display the general form obtained for the 
three-particle reduced density matrix of the NIFG, which 
will be the central quantity of our analysis. In Sec. Ill 
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we consider the classes of tripartite entanglement identi- 
fied for mixed three-qubit states by Aci'n et al. [l5j and 
discuss the properties of entanglement witnesses (EWs) 
capable of signaling the presence of these classes. A gen- 
eral scheme for constructing parametrized GHZ and W 
EWs via linear programming (LP) is introduced as a spe- 
cial case of convex optimization. Sec. IV is devoted to 
explicit development and application of new classes of 
parameterized EWs designed to detect genuine tripartite 
entanglement in the NIFG. First, adapting ideas from the 
work of Giihne et al. [l(| and Vertesi [l3j], we consider 
W entanglement witnesses based on a periodic spin-chain 
model. Second, we introduce a class of GHZ EWs which 
are constructed from projection operators corresponding 
to all the different types of state vector belonging to a 
system of three spin-1/2 fermions. Third, in order to 
identify the type of genuine multipartite entanglement 
that is present, we apply a stabilizer-operator formalism 
It is found that GHZ \ W genuine tripartite 
entanglement is not generated in the ground state of the 
NIFG as represented by the three-fermion reduced den- 
sity matrix, and that genuine W^-type tripartite entangle- 
ment, although present, does not exist in the absence of 
bipartite entangled and/or fully separable states. On the 
other hand, some genuine multipartite entangled states 
that have a positive partial transpose (PPT) with respect 
to some bipartition are found to occur. These conclusions 
are summarized in Sec. VI. Some details relating to the 
expectation values of the operators over the most general 
form of quantum states in the W set as well as the B set 
are collected in an appendix. 

II. INDISTINGUISHABILITY AND THE 
NONINTERACTING FERMI GAS 

Consider a system consisting of n parties {Mi}f_ v A 
/c-partite split is a partition of the system into k < n 
sets {Si}i =1 , each of which may be composed of several 
original parties. A given state p £ £>("Hdi <8 ••• ® Hd k ) as- 
sociated with some fc-partite split is called m-separable 
if a convex decomposition of it can be found such that, 
in each pure-state term, at most m parties are mutu- 
ally entangled, these not being entangled with any of 
the other n — m parties. In particular, a 1-separable (= 
separable) density matrix operator p <G 13(H) (belong- 
ing to the Hilbert space of bounded operators acting on 
H = Wdi <& ■•■ <£)Hd n ) is fully separable, being expressible 
as 

p scp = J2pM 1] )(^ } \® \af ] ){af ] \® ... ® \a^){a^\, 

i 

(i) 

with pi > and Y^iPi = 1- The system is called en- 
tangled when the corresponding density matrix operator 
is not separable. According to these definitions, separa- 
ble states necessarily form a convex set, since any con- 
vex combination of separable states is again separable — 
which is not the case for non-separable states. Beyond 



bipartite splittings, many different types of entanglement 
among the parties are possible, even for the case of dis- 
tinguishable particles. 

A schematic model involving two electrons located in 
a double- well was discussed in Refs. H, @ to illus- 
trate the consequences of indistinguishability for entan- 
glement. (See especially the related treatment of Ref. 0.) 
The qubit is modeled by the spin degree of freedom (with 
states | f) and | 1} ) , and there are two spatial wave func- 
tions labeled \4>) and |x), initially localized in the left and 
right potential well, respectively. For this bipartite sys- 
tem in a pure state, the authors considered the density 
operator p w = \w)(w\ E A(C 2K ®C 2K ). Denoting by f a 
and f\ the fermionic annihilation and creation operators 
for single-particle states a = 1, . . . 2K constituting an or- 
thonormal basis in C 2K , the ket \w) can be represented as 

\ w ) = J2 a ,b w abflfb\ Q )' witn the w °-b = ~ w ba defining an 
antisymmetric matrix. For any complex antisymmetric 
nx n matrix (w a b), there exists a unitary transformation 
U such that w' = UwU T has nonzero entries only in 2 x 2 
blocks along the diagonal i.e., 

w' = diag[Zi, ...,Z r ,Z ], Z. t = ( ^ ^ , (2) 

where Zi > for i = 1, r and Zq is the (n— 2r) x (n — 2r) 
null matrix. Each 2x2 block matrix Zi corresponds to 
an elementary Slater determinant. The matrix w' en- 
ables an expansion of the ket \w) in a basis of elemen- 
tary Slater determinants with a minimum number r of 
non- vanishing terms, r being termed the fermionic Slater 
rank of \w). A Slater rank of at least two is required 
for qualification as an entangled state. While this model 
is illuminating, its extension to more than two particles 
becomes very cumbersome, obscuring the nature of the 
correlations involved. 

An alternative description [9MTTI [l3j , more fruitful for 
our development, places the emphasis on reduced density 
matrices of a noninteracting gas of many identical spin- 
1/2 fermions (NIFG). By the Pauli exclusion principle, at 
most two such particles, with different spin values s, can 
occupy the same momentum state p = hk. The ground 
state of the system can be expressed as 

\^) = Y[bi( P )\o), (3) 

where [6j(p), &t(q)]+ = 5 st 5(p — q) and |0) is the vacuum 
state. 

Although the state of the system is written as prod- 
uct state, there are specific Pauli-cxclusion correlations 
between the constituent fermions arising from the com- 
mutation rules of the creation and annihilation operators 

and b. For the bipartite and tripartite configurations 
relevant to our investigation, these correlations are made 
explicit by deriving the two-body and three-body (two- 
fermion and three-fermion) reduced density matrices of 
the pure state . The two-fermion reduced density ma- 
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trix is given by 

Pss>- t W = (0ol4(r')4(r)^(r')^(r)|0o), (4) 

where ip\(r) [ip s (r)] creates [destroys] a particle with spin 
s at the location r. With the transformation 



(5) 



it is straightforward to obtain the two-body reduced den- 
sity matrix in the form 

Pss'-w = n 2 {S st S s , tl - 5 sn S sV f 2 (\v - r'Q], (6) 

where n = kp/ 6ir 2 is the density of particle of a given 
spin, and 

I 



is known as the Slater factor. (We focus here on the 
NIFG in three dimensions.) 

To evaluate the bipartite entanglement corresponding 
to the two-fermion reduced density matrix, one performs 
a partial transposition and determines the eigenvalues of 
the density matrix and its partial transpose Q. Entan- 
glement exists for two-fermion configurations such that 
f 2 > 1/2, i.e., for < |r — r'| < r e , where r e is 
cut-off radius for entanglement (rather than the clas- 
sical correlation). For the 3D NIFG, r e is determined 
by j 2 (kpr e ) = 1/2, where ji is the first-order spherical 
Bcsscl function. 

The same steps may be used to derive the reduced 
three- fermion density matrix of the NIFG as a function of 
particle locations r, r', r" and spins s, s', s" . Six possible 
arrangements give rise to six terms: 



p(s, s, s"; t, t', t") =(0 O |V4" ( r ")i>t' ir')il)\{r)il) s {r')ij} a ' (r)t/v (r")|0o) 

=n z {8 st 8 s >t'5 s »t» - fi2 f 'i3<5 ' s t$s't"S s >t" - fizf2?,5st>5s't5s"t" 

— 1 12 f '23^ st" S s 't'S s "t + J 12 j '13 j 23<5 st'fis't"$s"t + /l2/l3/23^st"^s't^s"t')- 

I 



(8) 



The three functions /12, /13, and /23 carry the respective 
arguments |r — r'|, |r — r"|, and |r' — r"|. 



The 9x9 matrix defined by Eq. (|8]) takes the form 



As 



/ r] 






















r)- 


P13+P23 
r 4 


-P23 
4 , 





-P13 

4 














-p23 
4 


~ 1 P12+P23 

' 4 





-P12 

4 

















T)- 


P12+P13 

r 4 





-P12 

4 


-P13 

4 








-pl3 

4 


— P12 

4 


T)- 


P13+P12 
^ 4 




















-Pl2 

4 


T)- 


P23+P12 

r 4 
















-P13 

4 







P23+P13 

r 4 





\o 




















v) 



r 



(9) 



where ?/ = (!- p i2 - P13 - P23V8 and 



Py = 



fii fijfikfjk 



"2 + /y + /ffc + /jfc fijfikfjk 



(10) 



Lunkes et aZ. [ll| have provided the following general 
expression for the n-body reduced density matrix of the 
NIFG, 



Pn 



a-E^i+J E wi*s><*< 



ln-2 ' 



(11) 



which is constructed from biscparable entangled density 



operators. Here, = 2~ 1 / 2 (|01) - |10)) is the maxi- 

mally entangled singlet state of the pair ij, and \fij\ < 1 
for all ij. Consequently \pij\ < 1 and < ?/ < j, 
so we have \p\2 + P13 + P23 | < 1 HH ■ This indicates 
that the density matrix of n noninteracting fermions can 
be written in terms of antisymmetric density matrices 
of fermionic pairs, albeit not in a convex combination. 
Considering the aforementioned entanglement condition 
for two-fermion configurations, entanglement would be 
present when all f 2 factors are greater than 1/2. Despite 
the explicit form (fTT|) of the n-body reduced density ma- 
trix of the NIFG as a combination of biseparable states, 
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the existence of genuine tripartite entanglement in this 
system was established in Ref. [TU 

The following sections will build upon this important 
result. We shall formulate new classes of three-qubit en- 
tanglement witnesses (EWs) and demonstrate that they 
can detect the corresponding fermionic tripartite density 
matrix 



(a) 



P3 =^I + a|*i2>(*i 2 l 
+ c|* 23 )(* 23 |® 
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(12) 



of the NIFG. (It is convenient to relabel the quantities 
Pij of Eq. dTO]) through a = p 12 , b = pi 3 , and c = p 2 3-) 
Obviously, p 3 cannot possess genuine tripartite entangle- 
ment for the case of simultaneous positive values of a, b, 
and c, due to the definition of biseparable states. 

It is instructive to sample the behavior of the coef- 
ficients a, b, and c of Eq. (jT2"j) for the NIFG. For one- 
dimensional (1-d) configurations specified by the distance 
x between fermions 1 and 2 and the distance r between 
1 and 3, Fig. [T] shows a plot of b with respect to a and 
c under varying kpr. Fig. [2] provides a similar plot for 
two-dimensional configurations in which fermions 1 and 
3 are separated by r and fcrmion 2 is constrained to move 
on a circle of radius r centered midway between 1 and 3. 
It is found that genuine tripartite entanglement cannot 
be generated in the region 4.5 < kpr < 5. 

Quite apart from consideration of the entanglement 
properties specific to the NIFG, the coefficients a, b, and 
c in Eq. (|12p must satisfy the inequalities 



1 -+ 1 -(a + b + c)± 1 -V~o 



b 2 + c 2 — ab — be — ac > 0, 



a 2 + b 2 + c 2 > ab 
?7>0, 



be- 



ac, 



(13) 



imposed by the restriction of the eigenvalues of any den- 
sity operator to positive-definite values. Referring to the 
definition of biseparable states, we note that by the def- 
inition of biseparable states, only negative values of the 
coefficients a, &, and c can give rise to genuine tripartite 
entanglement in p%. 




FIG. 1: For 1-d configurations as specified in (a), the coeffi- 
cients a and b are plotted versus c in (b) at different values of 
the quantity Ufx. Plot traces change from black to light gray 
as kpx increases from 0.1 to 5.0 in steps of 0.1. For fcjr > 4.5, 
three p functions attain positive values, placing the configu- 
ration outside the regime of genuine tripartite entanglement. 



one of them being compact, there exists a bounded func- 
tional (identified here as the witness operator) that serves 
to separate the two sets. For example, an EW can be 
defined for entanglement class e as an Hermitian opera- 
tor W such that (a) Tr(Wp s ) > for all fully separable 
states p s and (b) there exists at least one entangled state 
p e which can be detected by the condition Tr(Wp e ) < 0. 
It will be relevant to later developments that this defini- 
tion in itself cannot distinguish between different kinds 
of entanglement for more than bipartite systems. 

The objective of the present work is the elucidation of 
the tripartite entanglement properties of the NIFG. This 
objective is pursued within the framework of the compre- 
hensive analysis and classification of mixed three-qubit 
states provided by Acin et al. [TH|. This classification is 
a generalization of that for pure three-qubit states. To 
introduce the necessary definitions and results, we first 
recall that the most general pure three-qubit state takes 
the form 



III. ENTANGLEMENT WITNESSES 



Hghz) =A o |000) + Aie iy |100) + A 2 |101) 
+ A 3 |110) + A 4 |lll), 



(14) 



The existence of an entanglement witness (EW) for any 
type of entangled state follows from the Hahn-Banach 
theorem [2(|. In essence, this theorem establishes that if 
C\ and Ci are convex closed sets in a real Banach space, 



with Xi > 0, £\A? = 1, and 6 E [0,tt]. For the three- 
qubit system, there are two types of locally inequivalent 
entangled pure states. These are the GHZ type, defined 
by the form (|14|) with nontrivial parameters (especially, 
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FIG. 2: For 2-d configurations as specified in (a), the coef- 
ficients a and b are plotted versus c in (b) at different values 
of Icfx. Particle 2 is considered to move on a circle of fixed 
radius r, with X12 = r cos (6/2) and 123 = r sin (6/2). The be- 
havior of the coefficients is shown with traces changing from 
black to light gray as kpr increases from O.f to 5.0 in steps of 
0.1. 



A4 > 0), and the W type, represented generically by 

\ifar) = A o |000) + Ai|100) + A 2 |101) + A 3 |110). (15) 

The analysis of Ref. [IH identifies four distinct classes of 
mixed three-qubit states. The corresponding sets, de- 
noted by S, B, W, and GHZ, are all convex and com- 
pact, each being formed as a convex sum of appropriate 
projectors. Specifically, 5* states are mixtures of product 
vectors; B states are mixtures of product and bisepa- 
rable vectors; W states in turn are mixtures of prod- 
uct vectors, biseparable states, and W vectors (fT5")) ; and 
GHZ states are mixtures of all of the previous vector 
types as well as GHZ vectors (fPo)) . These four sets are 
invariant under local unitary and invertible non-unitary 
transformations. Evidently, they arc nested according to 
S dB <zW <Z GHZ. 

Our development will also involve non-convex subsets 



such as B \ S, W \ S, W \ B, GHZ \ S, GHZ \ B, and 
GHZ\W, which exclude fully separable states and hence 
contain only entangled states. A three-qubit state is said 
to possess genuine tripartite entanglement if it does not 
belong to the class of biseparable states B; such a state 
resides cither in GHZ \ W or W \ B. 

The classification of mixed 3-qubit states introduced 
by Acin et al. [l5| allows for the construction of entan- 
glement witnesses — namely GHZ EWs and W EWs — 
that are capable of detecting states with genuine tripar- 
tite entanglement. Thus, W EW will denote an opera- 
tor W w such that Tt(WwPb) > holds \/p B G B, but 
for which there exists a pghz\b G GHZ \ B such that 
Tr(yiV Pghz\b) < 0, thereby discriminating between 
the sets B and GHZ \ B. Similarly, a GHZ EW is pro- 
vided by an operator Wghz such that Tt(WghzPw) > 
is satisfied for any W state, but produces a negative ex- 
pectation value for some GHZ \ W state. Addition- 
ally, genuine tripartite entangled states, belonging to 
GHZ \ W and W \ B, can be identified by means of an 
EW operator designed to detect W\B entanglement — as 
will be demonstrated in the next section. 

A general scheme utilizing linear programming (LP) 
to arrive at EW operators for the detection of GHZ \ S 
entangled states was introduced in Refs. [U [U To con- 
firm the presence of genuine tripartite entanglement in 
the NIFG one needs W or GHZ EWs. Construction of 
parameterized GHZ EWs via the linear programming al- 
gorithm proceeds as follows. First, consider an Hermitian 
operator W in the form 



(16) 



possessing at least one negative eigenvalue, where the Pi's 
are Hermitian linear operators such that x < Tr(Pipw) < 
y for all and for any density operator pw G W. 

The parameters a; £ 1 are to be determined such that 
W qualifies as a GHZ EW. As pw varies over W states, 
pW _ Tr(Pipvy) maps W states into a convex region, 
since a linear functional maps a convex domain (here, 
the W class) to a convex region. Our principal task is to 
choose proper operators Pi so as to obtain an approxi- 
mating convex polyhedron surrounding the convex region 
spanned by the P^ (the so-called feasible region for the 
LP optimization). The term "approximating" refers to 
the fact that in general, not all points in the convex poly- 
hedron are produced by the expectation values of the 
Pi's over W states. By using this approximating convex 
polyhedron, the expectation value of the operator W over 
all W states is non-negative and, with respect to other 
states, admits at least one negative value. Accordingly, in 
seeking to determine a GHZ EWs of type (|16p . one needs 
to find the minimum expectation value of Y2% a iPi^ over 
the feasible region. In this way, the problem is reduced 
to optimization of the linear function aiP^ over the 
convex set provided by the approximating convex poly- 
hedron. 



6 



In characterizing the properties of the density operator 
for the NIFG, it is of interest to investigate the possi- 
bility of positive-partial transpose (PPT) entanglement, 
specifically the presence of a PPT entangled state with 
respect to some bipartition of a three-particle subsystem. 
The decomposability or non-decomposability of an EW 
may depend on which particles are involved. By defini- 
tion, an EW W is partially decomposable with respect to 
the i-th party iff there exist positive operators P, Q such 
that W = P + Q T *, where i € {1, 2, 3} and T< stands 
for partial transposition with respect to i-th party. Ac- 
cordingly, an EW is called partially non-decomposable 
EW with respect to a given party iff there exists at least 
one PPT entangled state associated with that party, i.e., 
p Ti > when Tr(Wp) < 0. The relevance of these def- 
initions will become clear when we encounter EWs that 
must explicitly bear the label of the party (or particle) 
with respect to which it is decomposable. For the pur- 
pose of identifying a PPT state with respect to a given 
party, one may derive the following conditions for the 
parameters in p 3 : A positive value for p Tl requires 

a 2 + b 2 + c 2 + ac + be - ab > 0, 

2(a 2 + b 2 +c 2 + ac + bc- ab) 1/2 - 1< -a-b + c<l, 

(17) 

the condition p T2 > is met only if 
a 2 + b 2 + c 2 - ac + bc + ab> 0, 

2(a 2 + b 2 +c 2 - ac + bc + ab) 1/2 - l<-a + b-c<l, 

(18) 

and p T3 > implies 
a 2 + b 2 + c 2 + ac- bc + ab> 0, 

2(a 2 + b 2 +c 2 +ac-bc + ab) 1/2 -l<a-b-c<l. 

(19) 

The presence of a PPT entangled density matrix with re- 
spect to a given party is signaled by satisfaction the cor- 
responding necessary condition among (|17I19[) together 
with its detection by the corresponding EW operator. 
Evidently, only a partially non-decomposable EW related 
to the selected party can detect a PPT entangled state 
with respect to that party. 

IV. EW OPERATORS FOR THE NIFG 

In accordance with the strategy proposed in Ref.H^, we 
now apply the LP method to determine proper EW op- 
erators for identification of chosen entanglement classes 
in the density matrix of the NIFG. We first consider EW 
construction based on an extended spin-chain model, and 
then build a new class of GHZ EWs from an appropri- 
ate set of density operators. Finally, we utilize stabilizer 
operators to develop EWs that discriminate between dif- 
ferent classes of genuine tripartite entanglement. 



A. Spin-chain model 

Guhne et al. [lj| were the first to introduce an EW op- 
erator (W EW) to detect genuine tripartite entanglement 
based on a macroscopic spin-chain model. The explicit 
form of this EW is 

yy(gen) = (1 + ^/g)! + p 12 + p 23; ( 20 ) 

where Pjj = a^' ■ a^' , or more rigorously 
P 12 ■ 

P 23 =<7< 2 > . ff (3) 

=iv>4M a) + i (1 M 2) 4 3) + i (1 M 2 M 3) ■ (2i) 

The superscripts identify the parties involved and a = 
(a x ,a y ,a z ). As required, the operator yy( gen ) has pos- 
itive expectation values with respect to all states be- 
longing to B and has a negative expectation value with 
respect to some genuine tripartite entangled state in 
GHZ \ B. 

Implementing the EW of Guhne et al. [Hj], Vertesi 
[l3T ] has shown that in a particular three-fermion con- 
figuration, there exists genuine tripartite entanglement 
in the NIFG. Our goal is to identify the type of genuine 
tripartite entanglement present in the NIFG. We develop 
EW operators suitable for this system and seek boundary 
conditions for separated fcrmions exploiting the witness 
operators. 

Starting from spin-chain model used in Ref. [l6] to con- 
struct the EW of Eq. ([2H]) , we consider the operator 

Wi = c I + P12 + P23 (22) 

and ask whether the constant Co can be chosen so that Wi 
becomes a GHZ EW. To answer this question, we should 
find the minimum expectation value of Pi 2 + P23 over all 
states in W set. The domain of W states is spanned by 
the explicit general form (fTS)) for \ipw), subject to any 
local SU(2) transformation. Upon expanding such the 
general W vector as \^w) — k=o 1 Aijk\ijk), one 
can evaluate the coefficient and thus the expecta- 
tion values py? = w\PijY& w) ■ (See the appendix for 
details.) It is then easily seen that the witness Wi can 
only detect W\B states, since the eigenvector associated 
with the minimum eigenvalue —4 of Pi 2 + P23 belongs to 
W domain, so Wi can never be a GHZ EW. 

Evaluating the expectation value of (|20p with respect 
to P3, we obtain 

Tr(W (gcn V 3 ) = 1 + V5-3(a + c), (23) 

which cannot be negative and maintain the condition a + 
b + c < 1 specific to the NIFG. 

To investigate the existence of genuine tripartite en- 
tanglement, we adopt the more general form 

Pz - U P3 U^ (24) 
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E-h-o 



0.01 0.02 



0.04 0.05 0.06 



FIG. 3: Plot of the expectation value T = Tr(W (gcn) p3) ver- 
sus k,FX at k,Fr = 0.1, in the 1-d coordinate scheme of Fig. [1] 
Here, yy( gcn ) is the entanglement witness f|20f) and p' 3 is the 
transform (|24|l of the three-fermion reduced density matrix 
P3 under an arbitrary local unitary operator. Physically, the 
distance of particle 2 from particle 1 is varied while the dis- 
tance between particles 1 and 3 is kept fixed. The trace T lies 
inside the gray region when kpx varies between and kpr 
and the parameters a and f3 range over all allowed values. 



for the density operator, where 



U 



(3* a 
-a* j3 



(g>3 



(25) 



convex polyhedron spanned by the , i.e., the expecta- 
tion values of the Py with respect to the B class of states 
(as considered in the appendix) . As a first conjecture de- 
limiting the eigenvalues of the P,j . we propose the set of 
inequalities 



"3 < -Pfl) fflj P%3 — 1' 



_pB 
r l2 

_pB 
r 12 

pB 
M2 



pB 
^13 
pB 
^13 

pB 
^13 



^23 < 3, 
P?3 < 3, 
P% < 3. 



(28) 



The polyhedron so described does not encompass the re- 
gion spanned by the . However, by parallel shifts of 
the boundaries in relations ()28[) . one can find a proper 
approximating convex polyhedron for reduction of the 
problem of finding a W EW to one of linear program- 
ming. This step is outlined in the appendix. The result- 
ing optimization problem reads: 



performs an arbitrary local unitary transformation of the 
density operator p% of the NIFG. (The same transforma- 
tion is applied for all three fermions, with |a| 2 + |/3| 2 = 1.) 
We then find 



Tr(W( gon V 3 ) = l + V5+(2-26-7a-7c)|/3 2 a 4 | 



9 (a + c)|a 2 /3 4 | - (a + c)|a| b - 3 (a + 



(26) 



Minimize 



subject to 



a 



ai2-Pf2 + ai3-fl3 + a 23-f2 B i 



3 < Pl2i Pl3i P23 — 1 




P?3 < 1 
pB _ pB < 1 
-"13 ■ r 23 — 1 

p?3 < 1 h 



Vs, 

v/8. 



(29) 



For one-dimensional configurations specified as in Fig. [T] 
Fig. [3] shows the value of Tr(W (ser V 3 ) as a function of 
kpx for fci?r = 0.1. It is seen that W\B entanglement 
is present over the range 0.01 < kpx < 0.09. 

Acknowledging the indistinguishability of the fermion 
constituents of the NIFG, we turn to a more general EW 
construction based on a parametrized operator that su- 
perposes all three of the spin products Pij : 



W (sp) = a I + a 12 P 12 + a 13 P 13 + a 23 P 



23- 



(27) 



For equal values of the real parameters a\ 2 = 013 = a>23> 
W (sp ), this ansatz reduces to W (sp) = a I 8 + a 12 (P 12 + 
Px3 + P23) which resembles a spin-chain model with peri- 
odic boundary condition. We now consider the possibility 
that can provide a W EW or even a GHZ EW by 

examining the constraints on the parameters ao, ai2, ai3, 
and a 2 3- 

In establishing W^p) as a W EW detecting genuine 
tripartite entanglement, the essential task is to find a 



With regard to maximum eigenvalues, the expectation 
values of 



P23 — 


Pl2- 


Pl3, 


Pl3- 


Pl2- 


P23, 


Pl2- 


-Pl3- 


' P23 



(30) 



reach 5. Therefore, solution of the LP problem (|2"9"f de- 
termines a region in which to form y\A sp ) as a witness 
operator. Solution proceeds by considering the intersec- 
tion of the constraints expressed in (|29[) and finding the 
vertices of the convex polyhedron. Upon solving the LP 
problem, one can obtain the values of the parameters ao, 
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fli2, ai3 and a23 satisfying 

a + (3 - \/8)oi2 - 3ai 3 + a 2 3 > 0, 

a - 3ai2 + (3 - VS)a 13 + a 23 > 0, 

a - 3ai2 - 3ai 3 + (—5 + Vo)a23 > 0, 

a + (3 - Vs)ai2 + a i3 - 3a 2 3 > 0, 

a - 3ai2 + ai3 + (3 - V8)a 2 3 > 0, 

a - 3ai2 + (-5 + v8)ai3 - 3a 2 3 > 0, 

ao + ai2 - 3ai3 + (3 - Vs)a 23 > 0, 

ao + ai2 + (3 - V8)ai3 - 3a 2 3 > 0, 

a + (-5 + V&)ai2 - 3ai3 - 3a 23 > 0, 

ao + ai2 + ai3 + «23 > 0, 

ao + ai2 + a i3 - 3a 23 > 0, 

ao - 3a i2 + 013 + a 2 3 > 0, 

ao + ai2 - 3ai 3 + a 23 > 0, 

ao - 3ai 2 — 3ai 3 — 3a 23 > 0, (31) 

provided that the expectation value of y\A sp ) is positive 
for all quantum states in the W set. Additionally, for 
W (sp) to qualify as a W EW, at least one of the eigen- 
values Ei among 

Ei = a + ai2 + a i3 + a 23 , 
E2,3 = ao — O12 — 013 — a 2 3 

± (af 2 + af 3 + a 23 - a i2 ai3 - ai 2 a 23 - ai 3 a 23 ) 1/2 , 

(32) 

must be negative, i.e., min(£^) < for i = 1,2,3, thus 
imposing on the parameters a; a condition additional to 
those of Eq. ([31]). 

In order to distinguish between different classes of gen- 
uine tripartite entanglement, one would like to find a set 
of parameters such that becomes a GHZ EW. Fol- 

lowing the same pattern as above, we seek a proper poly- 
hedron spanned by the P^ (see appendix, Eq. (A-5)). It 
is found that the maximum value of —P^ — + P23 , 
-P™ + P™ - P™, and - PW - P™ reaches 5, which 
is the maximum eigenvalue of the operators in Eqs. pop . 
Hence W^ sp ' can never be a GHZ EW for any choice of 
the parameters a^. It may be noted that use of the form 
yy( s p) m construction of an EW operator is compatible 
with the idea of developing it in terms of the projection 
operators |^r.)^r.| <gi I^ fc ) present in the expansion (TT2"]) 
of P3 , since one can write 

l*«X*«l®I (fc) = \<I-v®.aU)). (33) 

However, the more general (and more successful) imple- 
mentation of this idea pursued in Sec. IIVBI takes into 
account projectors referring to the set GHZ \ B. 

Returning to the task of detecting genuine tripartite 
entanglement in the three-particle reduced density ma- 
trix of the NIFG, we try W EWs corresponding to the 



ansatz 

°< p) = (1 + v/8) I + P12 + Pu - P23- (34) 
For the trace of the product of }V^ P ' with p 3 we obtain 

Tr( W^ p) p 3 ) = (1 + v/8)-3(a + 6-c), (35) 

which must be negative to signal the existence of W \ B 
entanglement in p 3 . For the NIFG the maximum value of 
a + b — c is 1. It can then be checked that the constraints 
imposed on the coefficients a, b, and c for the NIFG ex- 
clude the possibility of a negative expectation value of 
for the NIFG density matrix operator. We are 
led to conclude that symmetric spin-chain EWs are not 
capable of detecting genuine tripartite entanglement in 
NIFG. 

On the other hand, if we consider the general case of a 
density-matrix operator of the form (fT2]) subject to the 
basic constraints (flU)) . but not specializing to the NIFG, 

the expectation value of W^ P ' ) over p 3 can take on a 
negative value for some set of parameter values. It should 

be emphasized that entanglement detected by in 
such a general density matrix p 3 necessarily belongs to 
the W \ B class rather than the other class of genuine 
tripartite entanglement, i.e., GHZ \ W . We note that in 
the case of the transformed density matrix p 3 introduced 
in Eq. ([2"3]). the relation 

Tr(°W^ p) p 3 ) =(1 + VS) + (-a - 36 + c)\a\ 6 
+ (-3a-3fe + 3c)|^| 6 
+ (-8a-% + 6c)|a 4 £ 2 | 
+ (-9a-9& + 9c)|a 2 /3 4 | < (36) 

must be satisfied for W \ B identification. 

Having explored periodic spin-chain models for EW op- 
erators, we next describe another approach to the prob- 
lem of constructing EWs for the NIFG via LP, based on 
projection operators making up the density matrix p 3 in 
Eq. (0. 

B. GHZ EWs composed of projection operators 

In this section, we develop a new class of proper GHZ 
EWs with the aid of projection operators (pure-state 
density matrices) corresponding to the different entan- 
gled density operators arising in tripartite systems. If 
an EW is to detect tripartite entanglement associated 
with a given class of biscparables (1-23, 2-13, or 3-12), 
or with genuine entangled states (W \ B or GHZ \ W), 
it must necessarily be constructed from operators hav- 
ing non-vanishing expectation value with respect to the 
class that is specified. In particular, as demonstrated in 
the Sec. ELS (see e.g., Eq. ([55])). a GHZ EW cannot be 
successfully built from projection operators of the GHZ 
set if projectors from the GHZ \ W class are excluded. 
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With this preface, we introduce the operator 
^ghz = aoh 

+ «i(i*r 2 )(*r 2 i ® i + I*i3)(*r 3 i ® 1 + 1 ® i^x^i) 

+ a 2 1 Wi ) (Wi I + a 3 1 W 2 ) <W 2 1 + a 4 ^23) <*i2 3 1 . (37) 



To reduce the task of determining parameter values ai 
such that Wg H2 becomes a GHZ EW operator to an LP 
problem, we need to find a feasible region. Basing a first 
conjecture on the maximum eigenvalues of the Pj we find 
that the extremum points 



which contains projection operators for all the different 
vector types involved for the three-fermion subsystem, 

|W 1 ) = ^=(|001} + |010) + |100» ) 

|W 2 } = -^(|011) + |110) + |101» ) 



l*r 23 ) = ^|(|ooo>-|iii)), 



(38) 



and the a^ are real parameters. First, to ensure that the 
operator (|3T|) qualifies as a GHZ EW, we should impose 
positivity of the expectation value of Wq HZ with respect 
to W states. We take the generic W state vector \^w) 
in the form (fT5l) . For simplicity we define the operators 



Pi :=2 l^ijX^jl® 1 - 

<i,j> 

P 2 :=S\W 1 ){W 1 \, P 3 :=3\W 2 )(W 2 \ 
h :=2|*r 2 3)(*r 23 l. 



(39) 



and therewith their corresponding expectation values 
P^ (* = 1)2,3,4) with respect to VF-vectors in terms 
of the coefficients explicated in the appendix. One 
readily finds that the maximum eigenvalue is 3 for Pi, 
P 2 , and P3 and 2 for P4. A straightforward calculation 
shows that P 1 M/ , P 2 W , and P^ can reach their maximum 
possible value, but for P)± one finds a maximum overlap 



(Pi 



W P W nff 



2 ,, 3 ,P?) =(0,0, 0,0), (3, 0,0,0), (0,3, 0,0), 



(0,0, 3,0), (0,0, 0,3/2) 



(40) 



cannot be vertices for a feasible region of our LP prob- 
lem. As checked numerically, some points lying outside 
the convex polyhedron with vertices specified by (|40[) cor- 
respond to negative expectation values for Wq HZ . To 
compensate, we extend the proposed region by a parallel 
shift of the boundary hypcrplane and reduce the problem 
to LP as follows: 



Minimize 



subject to 



(Mi 



4 

E 

i=l 



a,P 



< Pf , P 2 W , P 3 W < 3 



< PF < 2, 



0<pw + pi 



2Pf <f. 



(41) 

Imposing the latter constraints, the operator Wq HZ can 
still have one or more negative expectation values since 
the range of the expectation value of Pi + P 2 + P3 + 2P 4 
is bounded between and 4. Now we have a feasible re- 
gion formed by the intersection of the 4-dimcnsional rect- 
angular parallelepiped domain defined by the extremum 
eigenvalues of the Pi's and the constraining hypcrplane 



p w + p w + p w + 2 pw < 15 / 4 _ Wc in us t ra te the situ- 



ation in Fig. [4j projecting onto the P™ dimension. 
The resulting convex polyhedron has vertices 



{p w ^ p w ^ p w ^ p w ) =(0; Q) Qj 0)) (3; Q) Qj 0)) (3; Q) Qj 3/8); (3) Qj 3/4j Q); (42) 

(3, 3/4, 0, 0), (0, 3, 0, 0), (0, 3, 0, 3/8), (0, 3, 3/4, 0), (3/4, 3, 0, 0), 
(0, 0, 3, 0), (0, 0, 3, 3/8), (0, 3/4, 3, 0), (3/4, 0, 3, 0), (0, 0, 0, 15/8). 

I 

After solving the LP problem posed in (|41[) . we arrive at a set of the a^, 



ao > 0, ao + 3ai > 0, ao + 3a 2 > 0, ao + 3a 3 > 0, ao + 15a4/8 > 0, 
ao + 3ai + 3a4/8 > 0, ao + 3a 2 + 304/8 > 0, ao + 3a3 + 304/8 > 0, 
ao + 3ai + 3a 2 / > 0, ao + 3ai + 3a3/4 > 0, ao + 3a 2 + 3a3/4 > 0, 

a + 3a 2 + 3o x /4 > 0, a + 3a 3 + 3ai/4 > 0, a + 3a 3 + 3a 2 /4 > 0, (43) 

I 

imposed by the positivity of the trace of of the EW oper- ator over all W states. These constraints, together with 
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FIG. 4: Projection of the feasible region of the LP problem 
on the P^ dimension, showing the intersection domain of the 
region P^ + P^ + 2Pf < 15/4 and the domain of possible 
expectation values of P2, P3, P4, with respect to a generic tri- 
partite quantum state within the rectangular parallelepiped. 



the existence of at least one negative eigenvalue among 

a , a + 3ai, a + 3a 2 , a + 3a 3 , a + 2a 4 , (44) 

guarantee that Wq HZ qualifies as a GHZ EW. Among 
the combinations f|44[) . only ao is excluded from negativ- 
ity, with the others remaining available for signaling a 
GHZ \ W state. 

We are now prepared to search for GHZ \ W genuine 
entanglement in the NIFG. As an optimal case of our 
GHZ EW construction, the inequality in Eq. (|4ip spec- 
ifying the last boundary hyperplane of the present LP 
problem yields the explicit GHZ EW operator 



°W^= T I 8 -2 Y, l*«X*tf|-3|Wi><Wi 

<i,j> 

-3|W 2 )(W 2 |-4|^ 23 )(* r23 |. 



(45) 



Adopting the witness °Wq HZ for our search, 



we calcu- 
late its expectation value for the three-fermion reduced 
density matrix p 3 of the NIFG, obtaining the simple ex- 
pression 



(46) 



which would have to be negative to confirm the presence 
of GHZ \ W entanglement. However, as seen in Fig. [51 
the term r\ does not reach negative values, ruling out 
this possibility. Moreover, calculation of the expectation 
value of the best-case witness operator Wq HZ with re- 
spect to the rotated thrcc-fcrmion reduced density matrix 



FIG. 5: Plots of the term 77 in Eq. Q46p versus 9 (in radians) 
at different values of kpr for the 2-d configuration shown in 
Fig. [2] (Plot traces change from black to light gray as kF is 
increased from 0.1 to 5.0 in steps of 0.1.) 



of Eq. 124]) yields 



15 

T 

\a 4 P 2 \ |« 2 /3 4 | 

4 4 
1 1 e 10 



(-33 + 9a -36 + 9c) 



L T- + I£ T- ( n + a + b + c) 

4 4 



a *3p3 p*3 a 3 



(1 — a — b — c), 



(47) 



which again fails to attain negative values. 

In the next subsection, GHZ stabilizer operators will 
be employed to formulate still another class of entangle- 
ment witnesses. 



C. Stabilizer EWs 

Previous work has demonstrated the utility and ro- 
bustness of EWs based on stabilization operators (TtI — 
EH, EH- By definition a stabilizer operator S for state \ip) 
has the property S\ip) = \ip). Toth and Giihne El 
have shown that if some of the stabilizing operators for 
a given state are available, entanglement conditions may 
be found that detect states in the neighborhood of this 
state. Here we apply the stabilizer formalism to obtain a 
parameterized EW for detecting quantum states close to 
a GHZ entangled state. Also employing stabilizer oper- 
ators, we are able to discriminate between different types 
of genuine tripartite entanglement that could be present 
in the three-particle reduced density matrix. 

We begin by considering a linear combination 



yy(stab) = bQl + &i £ + ^ + b 3 



(48) 



of GHZ stabilizer operators Si 



( 1 ) n .(2) T (3) 



_.(1)_.(2)_(3) a 
&x <Jx &x , 02 



cri 'a 



and Si 2 = Si x S2, where the bj's are real 
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parameters. First, to find ranges of the parameters bi 
such that y\/ stab ) becomes a W EW, we must determine 
the domain spanned by the expectation values of the 
corresponding operators Si over the biseparable set of 



states. Straightforward calculation shows that this task 
can be reduced to the solution of the following LP prob- 
lem (see the appendix): 



Minimize b a + &iSf + b 2 S 2 + &12SJ2 

subject to I ^ ilS h + l_^ 2S % + l_^ il+i2+ ^ 2 < 1; V(* 1; , 2 ) e {0, 1} . (49) 

I 



It is to be noted here that the boundaries associated with 
separable states cannot exceed unity (22j . The operators 
corresponding to the second set of constraints in (1491) . 
i.e., 

(-i) ll ^ 1 + (-i) i2 5 2 + (-i) il+i2+1 5 1 .2, (ii,*2)e {0,1}, 

(50) 



are positive operators and thus cannot form EW opera- 
tors. In solving the LP problem stated in (|4T))) . we find 
that the constraints 



&o + y/2(h +b 2 - 6i, 2 ) > 0, 
60 + \/2(-6 x + b 2 + 61,2) >0, 





60 + 61 


- 62 - 61,2 > 0, 




60 - 61 


+ 62 -61,2 > 0, 


&C)H 


- 726i + ^ 


2^2+- 1 +^6l,2>0, 

f V2&2 + ^61,2 > 0, 
f ^b 2 + V2b h2 > 0, 


60 H 


- - 1+ 2 ^h ■ 


&cH 


- -^1 ■ 


&0 


+ 1 -fb 1 4 


- ^b 2 - v/26i, 2 > 0, 

- V2b 2 + ±^&i,2 > 0, 


&0 


+ - 


60 


- + 1 


~fb 2 + x ~fb x . 2 > 0, 



60 + V2(bi -b 2 + bi. 2 ) > 0, 
b - V2(bi +b 2 + b h2 ) > 0, 
60 - b\ - b 2 + 6i >2 > 0, 
60 + h + b 2 + bi, 2 > 
b Q + V26i + ^^b 2 + ^61,2 > 0, 
bo + + V2b 2 + =^h, 2 > 0, ( 51 ) 

bo + =^bi + ^b 2 + V2b h2 > 0, 
b + + =±±V2fc 2 - s/2b 1>2 > 0, 

b a + =±£&b 1 V2b 2 + =±£&h,2 > 0, 
b - V2h + ^^b 2 + =^h, 2 > 



guarantee positivity of W^ stab ^ over all biseparable states. 
However, for yy( stab ) to qualify as the W EW we seek, it 
also must possess at least one negative eigenvalue from 
the possibilities 

&o+(-ir&i+(-ir 2 &2+(-ir i+ ' t2 &i, 2 , e {o, 1} 2 

(52) 

Next, to enable discrimination between different kinds 
of genuine entangled states, a GHZ EW is required. Ac- 



cordingly, we should find values of the b coefficients in 
Eq. (|4"5|) such that yy( stab ) is positive over all the states 
in the W class and yet has at least one negative eigen- 
value. To this end, we search for a polyhedron spanned 
by the expectation values , , and SY 2 , which are 
functions of the coefficients Ay/, in the general W vec- 
tor \^w) as expressed previously. Following the same 
pattern as before, we are led to the LP problem 



Minimize b + biS? + b 2 S? + b 12 S^ 2 

t. , , r (-i) i ^r + (-i) i2 ^ 2 M7 +(-i) ii+i2 ^r2 < 2.98 w ,. . . rn 112 

subject to | \-iyisY+ \-l) l2 Sf + (-lW^ +*s w < 1 ' ^ %2) E *- ' ' ( > 

I 

In this case, operators corresponding to the first cluster of the constraints in (|53p . i.e., 

(-l) ll 5i + (-l) l2 S 2 + (-l)' 11+42 5i, 2 (54) 
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can serve as GHZ EW operators. Solution of the LP parameters, namely 
problem ([5^]) yields a new set of constraints on the 6 



60 + 2.98(61 + 62-61,2) >0, 
6 + 2.98(-6i + 6 2 + 6i >2 ) >0, 
60 + 61 - 62 - 61,2 > 0, 
60 - 61 + 6 2 - 61.2 > 0, 
60 + 2.986i - 0.996 2 + 0.996i i2 > 0, 
6 + 0.996i + 2.986 2 - 0.996i' 2 > 0, 
6 + 0.996i - 0.996 2 + 2.986i.2 > 0, 
6 - 0.996i - 0.996 2 - 2.986^2 > 0, 
6 - 0.996i - 2.986 2 - 0.996 M > 0, 
6 - 2.986i - 0.996 2 - 0.996^2 > 0, 



60 + 2.98(61-62 + 61,2) >0, 
60-2.98(61 + 62 + 61,2) >0, 
60 - 61 - 6 2 + 61,2 > 0, 
6 + 61 + 6 2 + 61.2 > 
6 + 2.986i + 0.996 2 - 0.996i,2 > 0, 
6 - 0.996i + 2.986 2 + 0.996i' 2 > 0, 
6 - 0.996i + 0.996 2 + 2.986i.2 > 0, 
6 + 0.996i + 0.996 2 - 2.986i.2 > 0, 
6 + 0.996i - 2.986 2 + 0.996 M > 0, 
6 - 2.986i + 0.996 2 + 0.996^2 > 0, 



required for yVA stab ' to have only positive expectation 
values over W class. 

To construct an EW operator y\A stab ) which can detect 
a genuine entangled state belonging to W \ B but not 
accept a state from GHZ\W, the parametric constraints 
(l5"Tj) should be satisfied, while some constraint among the 
set (f55j) should be violated. To fulfill the requirement for 
detection, at least one of the eigenvalues in the set (|52|) 
must be negative. 

Among the above EWs of the form (|48p . we try the 
following 



w (stab) = ^2 I + 5i + 5 2 - 5i2, 

°W { ghz = 2-98 1 + ^+52-512. 



(56) 



To test for a genuine tripartite entanglement in the 
NIFG, we evaluate the expectation value of the first of 
these operators with respect to the reduced density ma- 
trix P3. For GHZ \ B entanglement to be present, the 
quantity 



Tr(°W^ tab) p 3 ) = v / 2-2 C 



(57) 



must be negative, i.e., a > l/s/2. Referring to Fig. [2j 
this would occur for negative 6 in the 2-d configuration. 

To detect GHZ \ W entanglement in the NIFG, the 
quantity 



Tr( W^ 



p 3 ) = 2.98- 2a 



(58) 



should reach negative values, but this is not possible since 
\a\ < 1. Considering the EWs defined in Eq. (|5o]) . we see 
that for a W \ B entangled state to be detected, the ex- 
pectation value of Si + S-2 — Si 2 must lie between y/2 
and 2.98; for the NIFG this requires that the condition 
— 1.49 < a < — v2/2 is fulfilled. Importantly, upon refer- 
ring to Fig. [6l we confirm that GHZ \ W entanglement 
does not exist in the three-fermion density matrix p 3 of 
the NIFG, and that all the states detectable in p% by the 
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FIG. 6: Plot of the trace T' = Tr{(Si + S 2 - 5i 2 )p 3 } versus 6 
(in radians) at different values of kpr for the 2-d configuration 
shown in Fig. [2] (Plot traces change from black to light gray 
as fcpr is increased from 0.2 to 5.0 in steps of 0.2.) 



entanglement witness "W^' belong to the GHZ \ W 
set. 

A measure called negativity has been employed in 
Ref. 11[ in studying the entanglement properties of the 
NIFG. This quantity is defined for a trio of fermions 
i,j,k as N [i>jk] = (Hpjplli - l)/2, where Hpiplli is the 
trace norm of the partial transpose of the three-fermion 
reduced density matrix p^ of fermion i with respect to 
fermions j, k. One of the advantages of working with 
EW operators rather than the negativity measure is in 
the identification of PPT entangled states. Using the 
condition (|19|) on the coefficients a, 6, and c in the NIFG 
expression for p^ together with the entanglement witness 



°W^ t , one can determine the configuration domain for 
which there exists PPT genuine W entanglement with 
respect to the third fermion. The behavior of the trace 
T" = Tr{°yV^ tab ^p 3 } at three values of kpr is shown 
in Fig. [7] for both 1-d and 2-d configurations. Hence 
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FIG. 7: Trace T" = Tr{°W^ tab) p3} for detection of PPT 
genuine tripartite W \ B entanglement with respect to the 
third particle, (a) The trace T" is plotted versus Ufx for the 
1-d configuration shown in Fig.[T]at fcpr = 3.0 and fcjrr = 3.5. 
(b) Trace T" is plotted versus for the 2-d configuration 
shown in Fig. U for 3.0 < k F r < 4.0 in steps of 0.25. (Plot 
traces change from black to light gray as kpr is increased from 
3.0 to 4.0.) 



Wm is a partially non-decomposable W EW for the 
third party since it can detect a PPT state with respect 
to the third fcrmion. 

It is instructive to note that the eigenvectors corre- 
sponding to the eigenvalues of p 3 in (fl"2"j) belong to the 
W\B class rather than GHZ\W, affirming the absence 
of genuine GHZ \ W entanglement in the NIFG. By way 
of proof, we first rewrite the three-body reduced density 
matrix in the standard form p% = J^. cJj|a!j}(a,|, in which 
the di are necessarily nonnegative. Then we assume there 
exists a GHZ EW operator Wghz such that 



Tv{W G hzP3) =Y J d i lk ( W GHz\a i ){a i \) 



(59) 



is negative. This assumption is contradictory if there is 
no GHZ \ W contribution to the eigenvectors of p 3 . 

One question still remains (cf. Ref. Q): can there ex- 
ist tripartite VF-type entanglement in the NIFG, without 
any admixture biseparable or separable components, i.e., 



a state belonging purely to the W\B subset? This aspect 
can be investigated by establishing an upper bound on 
the trace of the product of a general density operator p w 
possessing W\B entanglement, and the three- fermion re- 
duced density operator p% . By virtue of Fermi exchange 
antisymmetry, it is sufficient to work with the operator 



e wr =|Wi)(W 1 | 



=3(1001} 



|010) + |100))((001| + (010| 



(1001 



together with its transform 



Q 



IW 



(5* a 
-a* (3 



Q 



w 



f3 -a 
a* (3* 



(g>3 



(60) 



(61) 



under an arbitrary local unitary transformation, which 
again yields a purely W-type genuine entangled state. 
We have verified numerically that the value of Tr(g' w p^) 
is always less than one; hence the three-fermion reduced 
density operator of the NIFG cannot take the "pure" 



form g' 



This finding is confirmed for the explicit 



sxpl: 

forms of p(s,s',s";t,t',t") and ps given in Ref. and 
reproduced in Eqs. ([8]) and ([9]). Consequently, one can 
conclude that in the NIFG, genuine tripartite entan- 
glement only occurs in the three-fermion reduced den- 
sity matrix p^ in the company of 1-separable and/or 2- 
separable entanglement in other partitions: one cannot 
generate "pure" tripartite entanglement in the nonintcr- 
acting Fermi gas. 



V. CONCLUSION 

We have introduced new classes of entanglement wit- 
nesses (EWs) for the purpose of identifying genuine tri- 
partite entanglement, i.e., W \ B and GHZ \ W states, 
in the three-fermion density matrix of the noninteract- 
ing Fermi gas (NIFG). We have reduced the task of con- 
structing suitable EWs for this system to well-defined 
problems of linear programming. Considering EW op- 
erators inspired by a periodic spin chain model, EWs 
composed of the projection operators over the different 
classes of tripartite systems, and GHZ stabilizer opera- 
tors, we have found that the genuine tripartite entangle- 
ment present in the NIFG belongs to the W \ B class. 
This result is confirmed in the structure of the eigenvec- 
tors of a general three-fermion reduced density matrix p^ 
of the NIFG. We have seen that genuine tripartite en- 
tanglement does not occur in "pure" form, but instead 
it appears mixed with B or S components. Addition- 
ally, using a partially non-decomposable EW, we have 
been able to detect PPT genuine W \ B entanglement 
with respect to the third party of a fermion trio in the 
NIFG. The general approach followed in this work can 
be applied to investigate multipartite entanglement in 
other quantum many-particle systems, whether consist- 
ing of fermions or bosons, possibly with higher spins, and 
whether the particles are interacting or noninteracting. 
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Appendix 

To obtain the most general form for quantum states in 
the W set, we consider the explicit form of \ipw) given in 
Eq. (fl~5|) . together with an arbitrary local SU(2) trans- 
formation applied for the different parties i according to 
|0>< -> Oi\0)i + and -> (3*\0} t - a*\l) t , with 

| | 2 + \f3i\ 2 = 1. We next write the general locally trans- 
formed W vector as \^w) '■= V\tpw), where 



v=(R « ?W « ^V(a-i) 



Rewriting the general W vector as an expansion \^w) = 
Ylij k=o l Aijk\ijk), the coefficients A t j k arc determined 
as 

Atoo = A Q!iQ!2a3 + XiPla 2 a 3 + X 2 Pla 2 Pl + X 3 (3^/3 2 a 3 , 
Ami = X aia 2 f3 3 + X 1 f3la 2 f3 3 - X 2 (3la 2 a* 3 + X 3 /3l/3 2 p 3 , 
A 010 = \ a axp 2 a s + + X 2 (i{(i 2 (3* 3 - Ag^a^ag, 

Aq U = Aoax/32/33 + AiftAA - A 2 /3*/3 2 a* - X 3 fta* 2 /3 3 , 
-4ioo = A ^ia 2 a 3 - Xiala 2 a 3 - X 2 a\a 2 Pl - X 3 al/3 2 u 3 , 
A W i = X f3ia 2 (3 3 - Xiala 2 f3 3 + X 2 ala 2 a* 3 - X 3 a{P 2 P 3 , 
A110 = Xo(3i(3 2 a 3 - Xialf3 2 a 3 - X 2 alf3 2 f3 3 + X 3 a\a 2 a 3 , 
Am = Ao/3i/3 2 /?3 - Xial(3 2 (3 3 + A 2 a^ 2 «5 + X 3 a{a 2 f3 3 . 

(A-2) 



For the expectation values of the operators Pjj of 
Eq. (f2Tj) over the W set, we now have 

Pg :=^w\Pi2\^w) 

= |ylooo| 2 + |^ooi| 2 -|Aoio| 2 -|A 011 | 2 
-|A 10 o| 2 -|^ioi| 2 + |Ai 10 | 2 + |A m | 2 
+ (A i A* 00 + A 011 A* 01 + c.c), 

= |A oo| 2 -|^ooi| 2 -|Aoio| 2 + |A 011 | 2 
+ |A 10 o| 2 -|A 101 | 2 -|A 110 | 2 + |A m | 2 
+ (-4ooi^Sio + ^101^*10 + cc), 

= |A oo| 2 -|^ooi| 2 + |Aoio| 2 -|Aoii| 2 
-|Ai o| 2 + |A 10 i| 2 -|Aiio| 2 + |A m | 2 
+ (-^ooo-^ioi + AooiA 100 + ^oio^4*n + AonA* 10 
+ c.c), (A-3) 

while for the M^-set expectation values of the stabilization 
operators Si's appearing in Eq. (|48|) we obtain 

SY :=(Aooo^4in + AooiA* 10 + AoioA* 01 

+ A on A* 00 + c.c), 
Sr :=|A 000 | 2 + |A 001 | 2 - |A 010 | 2 - |A on | 2 

~\Ai Q0 \ 2 ~\Ai i\ 2 + \Aii \ 2 + \Aiii\ 2 , 
SY,2 := (~ AoooAJ 10 ~ -4001-4*11 + -4oio-4^ 00 

+ A01L4; 01 + A W0 A* QW + A W iA* Q11 

- A 1W A* 000 - AmA* 001 + cc). (A-4) 

Furthermore, for the operators Pi's entering Eq. (TTB| . the 
expectation values over the states of the W class read 

p x =2<*TH(i*r2><*ul® I+ l*r 3 )(*r 3 l ® i 

+ 1® I*23>(*2 3 I)I*W) 

={A\ m + A* ow + A* wo ) x (A 00 i + A QW + Aim), 

p 2 w =m w \wi) 1 2 

=(A* mi + A* ow + A* wo ) x (A 00 i + A QW + Aim), 
P™ =3|<* w |W a ) | 2 

= (A*ou + A* 1W + A* W1 ) x (Aon + A 1W + A 101 ), 

pr=2K*r 23 l^)l 2 

= (A* 00Q - A* 1U ) x {A Q00 - Am). (A-5) 

Similar relations are generated for the expectation val- 
ues of the operators P^ , Sf, and Pf. 

To find a domain of the parameters in Eq. ((27)) such 
that W (sp) qualifies MS i\ W EW, we search for a convex 
polyhedron which is embedded in the domain spanned 
by the . To illustrate how the problem of finding a 
W EW based on >V^ P ^ is reduced to the LP problem 
stated in (|2l?|) . one of the boundaries specified in Eq. (f2T))) 
for the convex polyhedron is determined as follows, the 
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pattern for the others being similar [16[. With \ipB) = 
I ^1)1^23) denoting an arbitrary biseparable state having 
entanglement among the second and third parties, we 
have 

\-p&-p& + p&\ = 

\{M<®<® + + + a (S) + ,(3) + ,(3) 

+ ^)43) +(r (2) £7 (3) +(7 (2) (r (3)|^ )| 

< + 4 2) + 4 2) + 4 3) + 4 3) + ^IV^I 

< I + VS. (A-6) 

The Schwartz inequality has been invoked in the sec- 
ond step. Parallel considerations apply in developing the 



other EWs introduced in this paper. In particular, for 
one of the boundaries associated with the stabilizer EW 
corresponding to Eq. (|49| we obtain 



I Si + S2 S12 1 — 

\(M - oMotoa® + ofPa® + oV>o&>aM\1, B )\ 

< \(M^ } 4 3) \^b)\ + \(M^ } \^b)\ 
+ \(^ B \a^a^\^ B )\ 

< y/2. (A-7) 
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